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Abstract 

Cathode emission current densities and anode 
work functions are necessary inputs to the analy- 
sis of thermionic converter systems. These in- 
puts must agree with existing experimental results 
and be free of ambiguity. A simple empirical 
equation giving the emission current density in 
terms of the cesium atom arrival rate and the sur- 
face temperature is provided here. The Swanson- 
Strayer correlation for the maximum change in work 
function due to cesium adsorption is discussed. 


least can be used as complex data fits, straight 
forward empiricism can yield a simpler formula- 
tion. 

We present here a simple and accurate 
converter surface physics formulation. The re- 
sult is an explicit equation for J in terms of 
the bare work function cp 0 , u, and T. Since the 
derivation is purely empirical, agreement with 
experiment is improved. 

The Cathode 


Introduct ion 


Converter analysis would be facilitated by a 
simple equation giving the electron emission cur- 
rent density J explicitly in terms of the elec- 
trode temperature T and the cesium particle 
arrival rate q. A need to computer simulate a 
system containing hundreds of converters makes 
such an expression even mere desirable. 

The results must of course agree with exist- 
ing experimental data. Also the expression must 
be useable for predictions of J in cases for 
which there are no experimental data (e.g,, for 
[i larger than in existing data and for untried 
electrode materials). 

Although considerable progress has beer made 
toward a first -principles understanding of the 
effect of alkali metal adsorption on the electron 
work function*, there is as yet no basic calcula- 
tion of J as a function of T and q. Jhere 
are a number of semienpirical corre lat ions' } 
however, which relate J to T and cesium gas 
properties (usually the Cs reservoir temperature, 
T r ) . By far the most tested and discussed of 
these are the classical atomic physics approach 
of Rasor and Warner** and the molecular chemistry 
formulation of Levine and Gyftopoulos . 


Although both works have in scone instances 
enjoyed agreement with experiment, some basic 
assumptions of the Rasor-Warner^- 0- and Levine- 
Gyftopoulos^ 2- ^ correlations have been brought 
into serious question. Steiner and Gyftopoulos^* 
have refined in part the molecular chemistry 
approach of Levine and Gyftopoulos. Notable dis- 
agreement with experiment has been found for both 
Levine (or Steiner) - Gyftopoulos and the Rasor- 
Warner eorrela+.ionfi.i5"20 while these works at 


A good deal of information has been obtained 
empirically about converter cathode surfaces . 

For example, R. Breitwieser has shown^l that his 
experimental data is such that plots of InJ 
versus InP , where P is the cesium pressure, 
are nearly linear for a given emitter temperature. 
Also, it has long been recognized that cp versus 
T/T r plots (for those experiment' for which 
U = Pcs/ V 2 jtmkT r ' and P Cs = vapor pressure of 
Cs at temperature T r ), are nearly independent 
of T r and linear in the range 2.3 5 $ S 3*3 eV 
for the refractory metals , where cp is the 
effective work function*. Some experimenters' 22 
results show a scatter of at most 0.1 eV in cp 
versus T/T r plots. These workers provide alter- 
nate methods of data plotting to reduce data 
scatter. 

Recently, using the "plasma anode" techni- 
que 20 , Houston + has shown that data on the metals 
W, Pd, Ru, Os, Rh, Ir, and Pt taken at three 
cesium reservoir temperatures, 4l4 °K, 434 °K, 
and 454 °K, "fall on a single curve" when plotted 
in the form cp versus T/Tj.. Also his data are 
nearly linear for 2.3 5 <P _ 3«3 eV. Thus, it is 
tempting to write 

cp = a T/T r + b (1) 

or 

J a 120 T 2 exp [-(a T/T r + b)/kT] (2) 

where a and b are to be determined empiri- 
cally. Eqs, (1) and ( 2) apply only to those 
cases where u = Pc s / V 2itmkT r ' . Generalizations 
will be discussed later. 


* cp = -kT In ( J/120T 2 ) 


* See Ref. 1 for a survey. 


+ A preliminary version was presented at the 
Thermionic Convertor Specialist Conference, 
San Diego, October 25-27, 1965* 
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If experimental data were available for all 
candidates for cathode surfaces, then it would be 
merely necessary to fit for a and b. But of 
course such data are not available, so it is 
desireable to also empirically write a and b 
in terms of cathode properties. The bare work 
function, ep 0 , is a property which distinguishes 
cathode surfaces , so one might try to find 
a(cp 0 ) and b(cp Q ). 

Since the work function of a surface is very 
sensitive to surface contamination, history of 
heat treatment, sputtering, bulk impurity content, 
etc., the following criteria should be satisfied 
by experiments used to empirically determine 
a(cp 0 ) and b(cp 0 ). First, it would be well if all 
samples were measured in the same tube. Secondly, 
the bare work functions snould be measured when 
the cesiated data are taken. Fortunately, there 
are now data available of this nature as taken by 
two experimenters, Houston 20 and Wilson^ - . 

There is enough disagreement between their re- 
sults, however, that one needs to make a choice 
between the two. Houston offers at least plausi- 
ble arguments for some inaccuracies in Wilson's 
results, so Houston's data are used here. Such a 
choice is not necessarily final, however, for as 
progress is made experimentally , the empirically 
determined values of a and b can of course be 
improved. 

A least mean square straight line fit was 
used to deterjnine a and b for the metals 
reported on by Houston (except for Pi, for which 
there was insufficient data), and the results 
are given in Table 1. The r.m. s. deviations of 
(a T/T r + b) from the data, 2*^ , and the magni- 
tudes of the maximum deviations, ^y, are seen 
to be quite small. 

Table 1 


Results of l.m. s. Fit for a and b 



a (eV) 

b( eV) 

^■ais( e V) 

Amnv(eV) 

range in cp 
compared 
(eV) 

Pt 

1.13 

-1.401 

0.029 

0.044 

2.10-3.41 

Ir 

1.33 

-1.87 

0.014 

0.024 

2.04-3.45 

Ds 

1.35 

-1.83 

0.014 

0.033 

2.02-3.59 

Rh 

1.42 

-2.09 

0.022 

0.052 

2.31-3.60 

Ru 

1.43 

-1.93 

0.031 

0.076 

2.06-3.55 

a 

1.48 

-1.88 

0.013 

0.026 

2.14-3.26 


* A *■ Experimental data point - (a T/T r + b) 


As intimated earlier, if one were to stop 
here, one would be limited to polycrystalline 
sur faces of the metals studied in Houston's ex- 
periment. In searching for a way then to re- 
late a and b to cp 0 , one might take note of 
Houston's plots of T/T r for a given cp versus 
<p 0 . As these plots are nearly linear, one could 
do a l.m.s. fit to obtain a and fj, where 

T/T r = a <p 0 + p (3) 


If this could be done for two values of cp, say 
and CP2 , then we would have that 

cp x = a (a x cp 0 + P x ) + b 

cp 2 = a (a 2 cp o + 3 2 ) + b (4) 

and the solutions for a and b as functions 
of cp Q would follow immediately. 

It turns out that rather good fits can be 
obtained using Eq. (3) for essentially all values 
of cp within the range ~ 2.20 eV to ~ 3.4 eV. 
Plots of the data including the l.m.s. linear 
fit for 3.3 and 2.3 eV are shown in figures 1 
and 2, respectively. In order to obtain more 
insight into the limitations of the method, the 
l.m.s. fit to the Pt group data was extrapolated 
down to cp 0 = 3.5 eV, for comparison with Hous- 
ton's 20 data on Th and Hf, taken in a 
separate experiment but using the "plasma-anode" 
technique. Although we may well not be interested 
in cathodes with bare work functions of less than 
say 4.3 eV, it is interesting that at cp = 3*3 eV 
the extrapolation agrees well with data on both 
Hf and Th, provided that extensions of the 
linear portions of the cp versus T/T r curves 
are used. At 2.3 eV, the agreement is good for 
Hf but not Th. 

If we let cp]_ = 3-3 eV and cpg = 2.3 eV, 
Equations (4) can be solved to yield 

- 1 
a = 0.176 cp Q - 0.158 

b = 3.3 - a [0.702 cp 0 + O.I93] (5) 

From the preceding, it appears that these results 
would be useful for 2.3 — cp 5 3-3 eV and perhaps 
for 4.3 5 cp 0 5 6.0 eV for Cs on metal substrates. 

As a test of the self-consistency of the 
above procedure, T/T r = agcpQ + £3 was fitted 
at cpg = 2.0 eV. Then afajcpg +“Pg) + b was 
compared with 2.8 eV, and found to agree to 
within ± 0.05 eV for 4.1 5 cp_ 5 6.5 eV. Addi- 
tionally, Table 2 gives the deviations L 
using Eqs. (5). They are larger than in Table 1, 
but the are all less than 0.1 eV. 

Table 2 


Accuracy Evaluation for Eqs. (5) 




Anax( eV ) 

range in cp 
compared (eV) 

Pt 

0.099 

0.17 

2.10-3.41 

Ir 

0.035 

0.058 

2.04-3.45 

Os 

0.016 

0.028 

2.02-3.59 

Rh 

O.089 

0.17 

2.31-3.60 

Ru 

0.054 

0.076 

2.06-3.55 

W 

0.021 

0.034 

2.14-3.26 


* A * Experimental data point - (a T/T p + b) 
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and 


But our work is not finished. Eqs. (1) and 
(2) as written are not useful for converter sys- 
tems design. First, for a given T, T r , and 
cathode surface, cp (and hence J) will generally 
still depend on the converter configuration 
(interelectrode .spacing, wall temperatures, etc.). 
Second, it would be well if the empiricism could 
be used in cases where the adsorbate is not com- 
ing solely from a vapor, e.g. , as in an impreg- 
nated cathode where adsorbate is supplied from 
the substrate. 

We can see how Eqs. (1) and (2) could be 
ambiguous (cp depends on more than T and T r ) 
by using the fact that cp and J are unique 
functions of T and u for a given cathode sur- 
face, where p is the cesium arrival rate. If, 
for example, we have an electrode immersed in a 
gas of temperature T g , pressure p„, and parti- 
cle mass m_, u = p / y 2 nrr g kT g . If Knuds en 
flow conditions apply, a = u 0 = evaporation rate 
at the cesium reservoir . But if continuum flow 
conditions apply, u = V Tr/T" u 0 . So for a 
given T and T r , different u's and hence dif- 
ferent J's are possible. 

Thus to generalize Eqs. (1) and (2) (and 
thus generalize the usefulness of our already 
found a's and b's), we need to write T r =T r (u) 
for the case of Houston's experiment. In Houston's 
experiment, PCs =P g and T g = T r . Approximat- 
ing Heimel's 2 ^ empirical expression for P Cs , 
which makes use of recent vapor pressure data, we 
have* 


J = 120 T 2 (M/C) 1,2856 a exp (-11605 b/T) (8) 

Eqs. (7) and (8), used in conjunction with Eqs. 
(5) and/or Table 1 will be useful regardless of 
whether the cesium arrives from a gas under say 
Knud sen or continuum flow conditions or even if 
it arrives from the substrate of an impregnated 
cathode . 

The Anode 

The temperature requirements on the anode 
(< 1200° K) are much less stringent than on the 
cathode. This, of course, gives us much more 
latitude. 

The anode generally has much more surface 
impurities than the cathode, so Houston's results 
are not necessarily pertinent to real convertor 
anodes. Still it would improve our under standing 
of real convertor anodes if one could predict the 
minimum work function of "clean" cesiated anodes. 

Although the minimum work function, 'p m ^ ni , 
is 1.5± 0.1 eV for many substrates , ^5 >2o it is 
desirable to have a means of predicting :%i n , 
generally, since we have some freedom in choosing 
anode materials. Swanson, 5 et al. have empiri- 
cally determined that 

AP m = 1.09 (1.78-<p 0 ) (eV) (9) 


M 3 C exp (-9027/T r ) 
so that 

1/T r = - (1/9027) In (p/C) (6) 

where C = I.I93 x 10 3 ° T r " 1,02 , 

T r = the mean T r in the range of interest , 

2 -1 

the units of p are particles cm sec , and 
T and T r are in °K. Values of C for use in 
converter design are given in Table 3, together 
with_the maximum error in u caused by the use 
of T r rather than T r in C. 

Table 3 

Arrival Rate Parameter C 


, -2 -1 
u( cm sec 

-2 -1 N 
C(cm sec ) 

Max. Disagreement 
with Heimel's u 

3.4xl0 17 -l.lxl0 19 

2.437x10 27 

8.956 

I.lxl0 iy -1.4xl0 20 

2.067X10 27 

7.6 % 

1. 4x10 20 -8 . 8xl0 20 

1.794x10 27 

6.6 % 


Combining Eqs. (1), (2), and (6), we have 
<P = -(aT/9027) In (m/C) + b (7) 


Other adequate vapor pressure relations are 
available (see, e.g. Ref. 26) 


or 

">0 + ^m = ®m = 1 '9 1+ - °-°9 <P Q (10) 

where cc^ is the minimum electron work function. 
Eq. (9) was determined by correlating field- 
emission data for Cs on Mo, W (Poly), W(100), 
W(110), Re, Ta, and Ni. Although thermionic 
measurements are done at much lower fields and 
higher temperatures, Table 4 shows that Houston's 
results are in good agreement with Eq. (9). Thus 
it appears that we may be able to use Eq. (10) 
for thermionic work functions also. 

Table 4 

A Test of Eq. (9) Under Thermionic Conditions 



<Po(eV) 

('Po-'Pmin. fteV) 

1.09(<p o -l-78)(eV) 

Pt 

5.50 

3.9 

4.05 

Ir 

5.31 

3.7 

3.85 

Rh 

5.25 

>3.7 

3.78 

Os 

5.16 

3.7 

3.68 

Ru 

4.89 

>3.3 

3-39 





W 

4.69 

>3-1 

3.17 

Hf 

3-97 

>2.3 

2.39 


* Ref. 20 
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Discussion 

It would be well to compare the accuracy of 
our scheme with that of the Levine (or Steiner )- 
Gyftopoulos and Rasor-Warner correlations using 
Houston's data. The Steiner-Levine-Gyftopoulos 
(S.L.G. ) schemes require knowledge of not only the 
bare work function, but also the monolayer work 
function and cesium coverage. As Houston did not 
measure these last two quantities, the S.L.G. cor- 
relations unfortunately cannot be compared unam- 
biguously. Houston 2 ^ compared the calculated 
(Levine-Gyftopoulos) S-ourve properties listed in 
Table 1 of Ref. 6 with his data (Table 1 was made 
up assuming a monolayer work function of„1.8l eV 
and a monolayer coverage of 4.8x10^ cm* for Cs 
on 21 metals). He concluded that "little corre- 
lation exists between the experimental results 
and the theoretical prediction." The predictions 
of the Rasor-Warner correlation are plotted on 
Figs. 1 and 2. There we see that at 3-3 eV, the 
Rasor-Warner predictions agree well with the data. 
At 2.3 eV, however, the Rasor-Warner correlation 
is in general disagreement with the data. 

Now let us see how our Eqs. (5) and (7) are 
used in the case of a convertor configuration. 

Ref eh et al. 2 ? have reported on a variable spac- 
ing convertor using a (110) oriented vapor de- 
posited tungsten emitter. They found the bare 
work function of the emitter to be 4.78 eV. The 
interelectrode spacing of 1 mil was such that at 
their highest cesium pressure Knudsgn flow condi- 
tions existed. Using C = 2.067xl0 2 ' cm" 2 sec _1 
(see Table 3) (In l/C = -62.90), and cpo = 4.78 eV, 
we can determine tp(u) from Eqs. (5) and ( 7 ). 

The result is shown to be in quite good agreement 
with the data in Fig. 3. 

It should be made explicit that the relations 
derived here do not include the effects of strong 
external fields. In taking this into account (as 
one might have to in ignited mode operation), it 
would be well to note Langmuir and Taylcr's 2 find- 
ing that the effect of external field on electron 
emission from cesium coated tungsten is in gen- 
eral larger than for pure tungsten and varies with 
coverage. 

Finally, in order to illustrate the use of 
the above equations, the zero-field electron 
emission characteristics are calculated for a 
convertor with the following properties: 

bare work function of the emitter, cp 0 = 4.9 eV 

bare work function of the collector, 

cpoc = 5.0 eV 

emitter temperature , T = I9OO ° K 

collector temperature, T„ - 850 °K 

* V 

u = 1.43X10 20 cof^sec - ^ at the emitter 

From Eqs. (5), a = 1.42 eV and b = -I.87 eV. 
Then from Eqs. (7) and (8) we have (using 
C * 1.794xl0 2 ^ cm“ 2 sec’l): 

9 * 3.02 eV and J = 4.35 amp cm -2 , 


where cp and J are of course the cathode 
effective work function and electron emission 
current density respectively. 

If we have a "clean" collector surface, 
Eq. (10) gives its minimum work function as 
<ftnin. = l**+9 eV - 

Summary 

An empirical convertor surface physics 
formulation has been provided which is quite 
simple and yet agrees better with experiment 
than other correlations in general use. 
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Figure 1. - T/T r versus substrate bare work function, 
<Pq, for (p ■ 3. 3 eV. 


<p, eV T/T r , (FOR <p = 2.3 eV) 



Figure 2. - T/T r versus substrate bare work function, 
<Pq, for <p « 2. 3 eV. 



Figure 3. - A comparison of the results of ref. 27 
with the predictions of eq. 7, using <p 0 » 4. 78 eV. 
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